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Abstract. In this paper a cellular automata model for one-lane traffic flow is 
presented. A new set of rules is proposed to better capture driver reactions to 
traffic that are intended to preserve safety on the highway. As a result, drivers 
behavior is derived from an analysis that determines the most appropriate action 
for a vehicle based on the distance from the vehicle ahead of it and the velocities of 
the two neighbor vehicles. The model preserves simplicity of CA rules and at the 
same time makes the results closer to real highway behavior. Simulation results 
exhibit the three states observed in real traffic flow: Free-flow states, synchronized 
states, and stop-and-go states. 



PACS numbers: 89.40.-a, 45.70.Vn, 64.60.My 05.40.-a 

1. Introduction 

Traffic networks are very complex systems where elaborated topologies are combined 
with large number of vehicles running on the network. Predicting traffic behavior is 
very important for planning and operation purposes. In the last years, computer 
simulations as means for evaluating control and management strategies in traffic 
systems have gained considerable importance because of the possibility of taking into 
account the dynamical aspects of traffic. 

In principle, traffic simulation models can roughly be divided into macroscopic 
and microscopic ones. While macroscopic models examine the dependencies between 
traffic flow, traffic volume, and average velocities; microscopic models investigate the 
movements of individual vehicles. In general, traffic flow models should keep the 
description of the relevant aspects of the flow dynamics as simply as possible by 
keeping track of the essential. In this spirit. Cellular Automata (CA) models for 
traffic flow were developed. Its main advantage is an efficient and fast performance 
when used in computer simulations, due to their rather low accuracy on a microscopic 
scale. These CA models for traffic flow are discrete in nature, in the sense that time 
advances with discrete steps, space is coarse-grained and properties of the CA can 
have only a finite, countable number of states allowing for high-speed simulations, 
especially when they are performed on a platform for parallel computation[Tl|2l[3]. 

The basic idea of CA models is not to describe a complex system with complex 
equations, but let the complexity emerge by interaction of simple individuals following 
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simple rules. Discrete space consists of a regular grid of cells, each one of which can 
be in one of a finite number k of possible states. All cells are updated in a parallel 
way, in discrete time-steps. The new state of a cell is determined by the actual state 
of the cell itself and its neighbor cells. The discrete nature of CA makes it possible 
to simulate large traffic networks using a microscopic model faster than real time. 
In 1992 Nagel and Schreckenberg proposed a stochastic cellular automaton model of 
vehicular traffic [4] , which was able to reproduce some empirically observed non-trivial 
traffic phenomena like spontaneous traffic jam formation!. This publication captured 
the interest of the physicists community and ever since there has been a continuous 
progress in the development of cellular automata models of vehicular traffic. Recently, 
extended CA models for traffic flow have been proposed to reproduce even more 
subtle effects, like free flow, spontaneous jam formation, synchronized traffic and 
meta-stability. These models incorporate anticipation effects, reduced acceleration 
capabilities and an enhanced interaction horizon for braking [5] El [7] . Due to their 
design, cellular automata models are very efficient in large-scale network simulations 

Eilllioilll]. 

In this paper, a new CA probabilistic model for traffic simulation that incorporate 
modifications to the NaSch model trying to reach a compromise between fast 
simulations and fidelity of results is presented. The key idea is to preserve simplicity 
of CA rules and at the same time make them closer to real drivers behavior. Rules 
of vehicle interaction in [4] are modified to better capture driver reactions to traffic 
that are intended to preserve safety on the highway. As a result, vehicle behavior is 
based on a safety analysis that determines the most appropriate action for a vehicle 
to take, based on the distance from the vehicle ahead and its velocity. This analysis 
is inspired by the work presented in [121 1131 114] for manual and automated highway 
systems. Although CA models can be applied to multiple lane highways, in this paper 
a one-lane highway with a ring topology is used. The goal is to show the ability of this 
CA model to capture the basic phenomena of traffic flow and preserve the simplicity 
and rapidity like the NaSch model. 

Simulation results presented conflrm that this CA model can reproduce most 
common regimes in traffic: free flow, synchronized flow and congested flow. The 
relations derived from the density/velocity and density/flow curves are in agreement 
with the empirical fundamental diagrams that describe these relations in traffic 
analysis. The influence of the variation of speed on the flow is also found to be a 
factor of great importance in traffic synchronization. 

The paper is organized as follows. A description of the proposed CA model is 
presented in Section 2. In Section 3, simulation results in the form of fundamental 
diagrams and other useful traffic representations are analyzed. Special emphasis on 
the newly introduced safety based rules is given. Finally, Section 4 contains concluding 
remarks and a summary of findings. 

2. Definition of the model 

The model presented here is a probabilistic cellular automaton. It consists of N 
vehicles moving in one direction on a one-dimensional lattice of L cells arranged in 
a ring topology. The number of vehicles is fixed. Each cell can either be empty, or 
occupied with a single vehicle that spans one or more consecutive cells. The velocity 

^The model hereafter is referred as NaSch 
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of a vehicle is constrained to an integer in the range {0, . . . ,v^ax}- In this paper, 
vehicles are allowed to occupy more than one cell. The speed limit, v„iax can be 
different depending on the kind of vehicle under consideration: trucks, cars, etc. For 
simplicity, in this paper only one type of vehicle is considered and therefore the same 
maximum velocity will be used for all vehicles. The integer velocity, that corresponds 
to one of the vehicle states in this CA, is related with the number of cells that a vehicle 
advances in one time step4 The other state, position, is related with the cell or cells 
that each vehicle is occupying. 

Typical length of a cell. Ax, in [1] is 7.5 m, where it is interpreted as the length of 
a vehicle plus the distance between vehicles in a jam. In this paper, three lengths for 
each cell of the lattice are analyzed: 5.00, 2.50 and 1.25 m, with the idea of showing 
the influence of this length on simulation results^ The time step (At) is always taken 
to be 1 s, therefore, transitions are from t t + 1. This time step is on the order of 
humans reaction time as pointed out in [15]. It can be easily modified. With these 
three values of Ax and At, v = 1 corresponds to a vehicle moving from one cell to the 
downstream neighbor cell, and translates into 18, 9 and 4.5 Km/h, respectively for 
cells of 5.00, 2.5 and 1.25 m. In addition, the maximum velocity is set to v^ax = 6, 
12 and 24 cells, respectively, which is equivalent to 108 Km/h in real-world units in 
all cases. 

On the other hand, to better capture driver reactions to traffic that are intended 
to preserve safety on the highway, our model includes three safety distances that a 
driver uses to decide to brake (brake distance ddec), accelerate (acceleration distance 
dacce) or keep his/her velocity (distance to keep dkeep)- Drivers' decisions are based 
on preserving safety. Vehicles velocity is not solely determined based on the distance 
to the corresponding vehicle in front, also considers the speed and the deceleration of 
the front vehicle. In this way, the cell lengths considered (grid granularity) in the CA 
are the unit of distance and will have an influence in the way normal maneuvers of 
deceleration are performed. 

For the coarser grid, that corresponds to a cell length of 5.00 m it will be assumed 
that the minimum allowable deceleration will be reached in one time step (M = 1). 
For a intermediate granularity, with a value of cell length of 2.50 m, this will occur in 
two times steps (M = 2). Finally for the finer grid, with cell length of 1.25, reaching 
the minimum allowable deceleration will take four time steps (M = 4)^. Emergency 
braking in all cases will have a value of —5.00 m/s^ and will be reached in one time 
step. Maximum acceleration will be 5.00 m/s^, 2.50 m/s^ and 1.25 mj , for cell 
lengths of 5.00, 2.50 and 1.25 m, respectively (corresponding to accelerating from 
— 100 km/h in 6, 12, and 22s, respectively). 

Due to the discrete nature of space and time in CA models, units of distance, 
velocity and time are normalized with respect to the length of each cell, Ax, and the 
time step. At. Therefore, units in position x denote the number of cell in the lattice; 
in velocity v, number of cells per unit time, and in time t, number of time steps4 

The nth vehicle is characterized by its position a;„(t) and velocity ?;„(t) at time 
t. Vehicles are numbered in the driving direction, i.e. vehicle n + 1 precedes vehicle 
n. The space gap (where s is the size of nth vehicle, expressed in number of cells) 
between consecutive vehicles is denoted by (i„(<) = a;„+i(t) — a;„(t) — s; for which it is 

^Provided that there are no vehicles that obstruct its forward movement. 

"^These quantities corresponds roughly to full, half and quarter vehicle length, respectively. 

^Note that M * cell length = 5.00 m. in all cases 

^For this reason, for example ii < d is used instead of f < d/A.t, because At = 1. 
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assumed that a vehicle position denote the cell that contains its rear bumper For the 
nth vehicle its safe distances ddea dacce f^nd dfceep are defined in the following form: 

M 

dacce„(t) = Y[[K(i) + l) Div M] + l][K(i) + l) Div M] 

+ [{vn{t) + l) Mod M][[K(t) + l) Div M] + \] (1) 
M 

~^[[{vn+i{t)-M) Div + Div M] 

~[{vn+i{t) - M) Mod M][[{Vr,+l{t) - M) Div M] + l] 

M 

dkeep^{t) ^ —[[Vn{t) Div M] + l\[Vn{t) Div A/] 

+ Mod M][K(t) Div M] + l] (2) 

M 

-^[[{vn+i{t)-M) Div A/] + l][K+i(i)-A/) Div M] 
-[K+i(t)-M) Mod A/][[K+i(t)-M) Div M] + l] 



t^decjt)- Y[[K(t)-l) Div + Div M] 

+ Mod Div A/] + l] (3) 

~^[[{vn+i{t) ~ M) Div M] + l][K+i(i)-A/) Div M] 
-[{vn+i{t) - M) Mod Af][[(w„+i(t) - Af) Div A/] + l] 

Expressions in equations ([T]) to ([3]) represent the safe distances the n-th vehicle 
must have with respect to its preceding vehicle if it is going to accelerate, (H]), keep 
its velocity, (U), or decelerate, ^ in the current time step. Operations (a Div h) 
and (a Mod h) denote the quotient and the remainder, respectively, resulting from 
dividing a by h. The basis to calculate these safe distances is to assume that the worst 
possible scenario after any of the these three basic maneuvers is performed corresponds 
to the vehicle in front applying full brakes |16j . 

There are two groups of terms in the right-hand side (RHS) of each one of the 
expressions given in equations ([T]) to ([3]). The first two terms represent the traveled 
distance by the vehicle n assuming that in the next time-step it finds that the vehicle 
n -|- 1, in front of it, has slammed the brakes; forcing vehicle n to also hit the brakes. 
The last two terms in the RHS of the expressions given in ([T]) to ([3]) are the traveled 
distance by vehicle n -f- 1 if at the current time-step it slams the brakes. 

To determine state transitions the following set of rules, which are applied 
simultaneously to all vehicles, is defined: 

RO: Calculate ddecdt), dacdt), and dkeepit) 
Rl: Acceleration. 

If dn(t) > dacc„{t), the velocity of the car n is increased by one, i.e., 
Vn{t + 1)^ min{vn{t) + 1 

''Vehicles space gap represents effective distance, corresponding to the number of empty cells 
between vehicles 
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R2a: Cruising. 

If dacc„ (t) > d„ (t) > dkeepn {t) , vclocity of vehicle n is kept equal with probability 
1- R, i.e., 

Vn{t + 1) ^ Vn{t) with probability 1 — R. 
R2b: Random braking. 

If rfacc„(i) > dn{t) > dkeepn {t) and Vn{t) > 0, vclocity of vehicle n is reduced by 
one with probability R: 

Vn{t + 1) ^ max{vn{t) — 1,0) with probability R. 

R3: Braking. 

If dkeep„{t) > dn{t) > ddeccn{t) and Vn{t) > 0, vclocity of vehicle n is reduced by 
one: 

Vnit) max{vn{t) - 1,0) 
R4: Emergency braking. 

If Vn{t) > and dn(t) < ddecc„{t), velocity of vehicle n is reduced by M, provided 
it does not go below zero: 

Vn{t + 1) — > max{vn — M, 0) 
R5: Vehicle movement. 

Each vehicle is moved forward according to its new velocity determined in rules 
1-4: 

{t + 1) ^ Xn{t) + V„{t + 1) 

Rules Rl to RA are designed to update velocity of vehicles; rule i?5 updates position. 
According to this, state updating is divided into two stages, first velocity, second 
position. The rationale behind rules Rl to RA is as follows. 

Rl: This rule postulates that all the drivers strive to reach the maximum velocity 
whenever possible. This is in agreement with other velocity policies, as it is the 
case with the greedy policy in [17j . 

R2a: This rule reflects the fact that drivers will try to keep their velocity if they 
perceive the distance with the vehicle in front as safe. 

R2b: This rules is introduced to model traffic disturbances that cause drivers to 
reduce their speed for no apparent reason. These can happen, for example, 
due to incidents along the highway that distract drivers. This random braking 
contributes to creation of traffic jams. 

i?3: This rules requires the driver to apply moderate braking when the spacing that 
separates his/her vehicle to the vehicle in front is becoming small. 

i?4: This rule stresses the approach taken in this paper: the most important drivers' 
decisions are related to safety. Thus, when according to its speed and the speed 
of the vehicle in front, the driver perceives an unsafe spacing, he/she will slam 
the brakes. When conditions for this rule are met, driver is in an unsafe situation 
that could lead to a collision if the driver in front slam the brakes. For this reason 
the driver of the vehicle will slam the brakes. If initial distributions of relative 
distances to the velocity are selected to satisfy at least ddec„, then vehicles will 
never have a relative distance such that this rule is activated. This rule was 
introduced to proof this fact and to allow also perturbations in the other rules 
trying to investigate their effect. 

There are several modification to the NaSch model. First and most important, it 
should be noted that velocity setting rules depend not only on the relative velocity of 
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neighbor vehicles, they now take into account tlicir relative distance. This modification 
was included to incorporate normal drivers behavior that base their driving decisions 
on both relative velocity and relative distance. It should be noted that safe distances 
given in equations ((1]) to ([3]) grow faster than linear with relative velocity of vehicles. 
This is in accordance with normal drivers spacing policies, for example, those based 
on vehicle following, constant time headway, etc. Acceleration and deceleration 
magnitudes for basic maneuvers are much smaller than those allowed in the NaSch- 
model. The values that are used here are much closer to those reported for normal 
driving. Emergency braking deceleration of —5m/s^ is considered an acceptable for 
this maneuver [14]. Note that the order for applying rules Rl to i?4 is not relevant as 
conditions for each one form disjoint sets. 

Other relevant modification to the NaSch-model is the change in the application 
of the deceleration and randomization rules. In the NaSch-model, randomization is 
applied after deceleration, while in the model here proposed randomization is applied 
only to vehicles that are in cruising conditions and do not require to brake. In this 
way double braking is avoided. The only probabilistic behavior is included in rule 
R2b. It should be noted, however, that the value of the probability of random braking 
will be smaller than that reported in [4]. This is consistent with the idea of treating 
random braking as a disturbance that should not occur very often. 

This CA model is a minimal model in the sense that all the five steps are necessary 
to reproduce the basic features of real traffic. However, additional rules may be 
incorporated to capture more complex situations [5]. The parameters of the model 
are the following: number of cells L, maximum number of vehicles N^ax, number of 
vehicles driving N , limit speed Vmax, vehicle length s, number of time steps to achieve 
maximum braking AI and the random braking probability in cruising R. 

3. Simulation results 

In this section, the experimental setup of the proposed model is introduced and results 
obtained are described. 

3.1. Simulations setup 

To simulate the CA model proposed in the previous section, the total number of cells 
is assumed to be L = N^^x * s, where N^ax denote the maximum number of vehicles 
that can be driving on a circular lattice and it is set to 10**. Each simulation starts 
with a initial configuration of N vehicles, with random distributions of speeds and 
positions. Since the system is closed, the density, p = N/L, remains constant with 
time. Based on this considerations, the maximum density will be 200 veh/km in 
real- world units. 

In order to prevent traffic accidents as a previous step, speed values are adjusted 
in such a way that dn > ddecc- Starting from this initial configuration, the average 
density, velocity and fiow, denoted by p, u, and q, respectively, over all cars are 
measured each time step. 

All the simulation data presented in this work was calculated using T ~ l5*Nmax 
time-steps. In order to analyze results, the first 10 *N„iax time-steps of the simulation 
are discarded to let transients to die out and the system to reach its steady state. 
Then the simulation data are averaged over the final 5 * N^ax time-steps. Traffic is 
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also considered to be homogeneous, so all vehicle characteristics are assumed to be 
the same. 

Velocities are updated according to the velocity updating rules R1 — R2 — R3 — R4: 
and then all cars are moved forward in step R5. For each simulation a value 
for parameter R is set and global densities, flows and speeds are calculated. In 
addition, a local detector with 200 * s cells is used to compute measurements of these 
same variables. In this latter case, data points were collected with a measurement 
period equal to 300 time-steps. Fundamental diagrams are constructed based on this 
considerations. 

3.2. Different length of cells and random braking level 

As mentioned above, in the proposed model three lengths for each cell of the lattice 
are considered (Ax). Determining the influence of this length on simulations results 
is important in order to propose parameters of the model that reproduce behaviors 
closer to those observed in the reality. Following that proposal, traffic flow behavior 
using the model of previous section with values of Aa; equal to 5.0, 2.5 and 1.25 m was 
investigated. Each of these values is associated to a different value of the parameter 
M (which represents the number of time steps to achieve maximum braking) ,1,2 and 
4, respectively. 



Ax=5.00 m 




50 100 150 200 

density(veh/km) 



Figure 1. Fundamental diagram for different values of cell length Ax. The 
smaller length cell is, the higher flow is. 

In figure [T] the fundamental diagram of the proposed model with a fixed value 
of i? = 0.15 and different values of Aa; is shown. This diagram characterizes the 
dependence of the vehicles flow on density. In this figure, the curve shows points 
for p varying from 2 to 200 veh/km in steps of 2 veh/km. From this diagram, the 
influence of the length Ax on accelerations and decelerations can be observed. Smaller 
values of Ax, that is lower acceleration levels, imply larger flows. This behavior can 
be explained in the following way. The emergency braking follows the same pattern 
independently of the values of Ax (the same maximum deceleration). Nevertheless, 
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with a finer grid the safety distance to accelerate required by a vehicle is smaller, 
due to the fact that the acceleration capacity is lower. This produces larger velocities 
leading to an increase of the flow. 



2500 -, 
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Figure 2. Fundamental diagram for cell length Ax = 0.15 and different values 
of the parameter R. 

It is interesting to note that the initial positive slope, corresponding to a free-flow 
state where there are no slow vehicles, is similar for all values of Aa; as maximum 
speed is independent of cell length. Here all vehicles travel with the maximum speed 
Vmax- On the other side of the plot, for all cell lengths, curves reach the same curve 
for large densities where the flow decreases with increasing density. From flgure [1] 
it is clear that the maximum flow for Aa; = 2.5m is closer to values reported in the 
literature for one lane highway [SI [6]. So, this value of Ax is adopted for all data 
presented in the rest of this paper. 

The next step was to investigate the influence of the probability of random braking 
R on the simulation results. Figure [2] shows the average fundamental diagram for 
different values of the parameter R. As can be noted from this flgure, for very small 
values of R there are oscillations in the fundamental diagram that are induced by the 
integer arithmetic. When R reaches a value of 0.10 these oscillations almost disappear. 
Note that for values of R from 0.10 to 0.20 there is almost no change observed in the 
fundamental diagram. Thus, value of R was set to 0.15 for the rest of the simulations 
presented in this section as with this value there are very small oscillations in the 
fundamental diagram and the corresponding value of its maximum flow is closer to 
that reported in the literature for a one lane highway. 

3.3. Traffic flow organization 

The proposed model is able to reproduce the three states observed in real traffic fiow: 
(i) Free-flow states, which are characterized by a large mean velocity, (ii) synchronized 
states, where the mean velocity is considerably reduced compared to that of free flow 
states, but all cars are moving, and (iii) stop-and-go states, where small jams are 
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Figure 3. Fundamental diagram resulting from random initial conditions. 
The simulations are performed on a ring with a length of 20000 cells, each 
corresponding to 2.5 m. The parameters of the model are Vmax = 108 km = 
12 cells s-^, R = 0.15, Ax = 2.5 m, and M = 2 

present. Synchronized traffic and stop-and-go states will be considered as congested 
states. 

In figure [31 a fundamental diagram obtained from random initial conditions for 
parameters values of i? — 0.15 and Ax = 2.5 m is shown. Each dot in the plot 
represents a measurement taken by the local flow detector described in the previous 
section. The white lane is the average flow during all the simulation. One can 
distinguish, in the fundamental diagram, three different states depending on the global 
density, namely low density range, high density range and intermediate range for 
pi < p < Ph. By increasing p in this intermediate range, the flow presents a slow 
decrease. The explanation for this behavior is related to the fact that according to 
rule RA of the proposed model drivers can decelerate in more than one unit. In this 
intermediate range, increasing the number of vehicles leads to relative distances that 
imply deceleration ddec^ to be smaller than the average space gap between two vehicles, 
^ — s. This forces rule RA to be applied (emergency braking). Vehicles that perform 
emergency braking produce an instantaneous block of the system and induce other 
vehicles to reduced their velocities. This causes a local decrease in flow. 

The shock waves and jams formations for these three flow regimes can be observed 
in time-space plots shown in figurelH The time-space diagram is a graph that describes 
the relationship between the location of vehicles in a traffic stream and the time as 
the vehicles progress along the highway. In this flgure, each vertical rows of dots 
represents the instantaneous positions of the vehicles moving towards the right of 
highway (top) , while horizontal rows are vehicles crossing the same highway positions 
at different time-steps. Trajectories of individual vehicles move forward time and 
space. One can observe the shock waves and jams formations. For the low density 
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Figure 4. Time-space diagram for R = 0.15, Ax = 2.5 m and different values of 
density. The solid lines in diagrams correspond to the trajectory for a specific car 
at different time-steps, (a) corresponds to a low density range, for p = 16 veh/km. 
(b)-(c) correspond to a intermediate density range for p = 30 and p = 50 veh/km, 
respectively, (d) corresponds to p = 140 veh/km 



range, these formations are moving forward (see figure |4lja)). After this range, the 
shock waves move backward (see figures|31Jb)-(d)), indicating the presence of congested 
states. Note that stopped cars with zero velocity only are presented in high density 
ranges (see figure Sl^d)). 

As it can be observed in velocity distribution shown in figure O for the 
intermediate density range, vehicles' velocity is considerably reduced when compared 
to the low density range, although all vehicles are moving, the situation represents 
a synchronized state. This figure also allows for a clear distinction between the free 
fiow and the congested state. In the former the mean velocity remains constant at a 
high value (vmax), as density increases, a capacity drop occurs resulting in a steady 
declination of fiow (thick solid line). Note that at the critical density, the standard 
deviation (thin solid line) jumps steeply; this means that vehicles' velocities start 
fiuctuating after the transition point. Once a compact jam is formed (after p^), the 
dominating speed quickly becomes zero. 

Figure [S] shows the results of a simulation with the same conditions, R = 0.15 and 
Ax = 2.5 m, as previous figures. In this case the emphasis is in platoon organization. 
Ordinates in figure [6] show the proportion of vehicles that are traveling in Vp-platoons. 
Here t;p-platoons are defined as a group of two or more neighbor vehicles traveling at 
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density[veh/krn] 



Figure 5. The distribution of vehicles' speeds, as a function of the global density 
resulting of the proposed model with _R = 0.15 and Ax = 2.5 m. The thick solid 
line denoted the space mean speed, whereas the thin solid line shows its standard 
deviation. The gray regions denote the probability density. 




density [veh/km] 



Figure 6. Mean number of vehicles by platoon with velocity Vp respect to the 
mean density 



the same speed, Vp, with the space between them dictated by rules Rl-RA. Figure [6] 
shows that for low densities all vehicles travel at maximum speed and therefore a single 
platoon is formed. For intermediate densities, there is a reduction in the number of 
vehicles that are traveling in platoons, although the total number is still significative, 
as about 40% of vehicles still travel in platoons. When density approaches the stop- 
and-go values, there is a marked drop in the number of vehicles traveling in platoons. 
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due to an increase in the dispersion of velocities (see figure [5]) • Finally, for very large 
densities, the number of vehicles in platoons increases as more and more vehicles 
have zero velocity. The gray levels in figure [5] are related with an specific velocity, 
according to the scale showed in its right hand side. It is interesting to note that 
vehicles traveling as platoons have very similar velocities. 

On the other hand, the proposed model in this paper is dependent on the initial 
conditions. In the next subsection this behavior is analyzed. 

3.4- Effects of the initial conditions 

In figure [7] two curves for the fundamental diagram for i? = in the range of densities 
of the interest are shown. The upper curve is calculated by starting from an a 
homogeneous state, where vehicles are distributed equidistantly with the same initial 
gap, {gini), and velocity (wmi); such that gini = ^ and Vini is equal to the maximal 
speed such that gi„i > dkeep- As opposed to the previous case, to obtain the lower 
branch the system was initialized by distributing the vehicles with random velocities 
and positions over the lattice and allowing the system to relax as time pass. 




50 100 150 

donsity[voh/km] 



Figure 7. Fundamental diagram resulting from two different initial conditions 
for R = 0. The upper curve was calculated by starting from an homogeneous 
state, whereas the lower curve was obtained from random initial distribution 

If system is initialized randomly, this has the effect that some vehicles are spaced 
more closely to each other. As a direct consequence of this, all vehicles should adjust 
their velocities according to rules i?3 — i?4. In the extreme case, a single vehicle 
with = 1 cell/s will cause that other vehicles behind apply emergency braking, 
resulting in total flow decrease shown in the lower curve in figure [71 When, the 
system is initialized homogeneously with the same gap and velocity for all vehicles, 
then between the densities pi and p2 a metastable high- flow branch exists. This is 
due to the fact that for these range of densities the average flow still depends on the 
initial configuration. For densities just below pmax the initial gap, gapini, allows to 
keep the maximum initial velocity equal to Vmax {go-Pmi > diueep, where diueep is the 
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minimum safety distance that should exist between two vehicles to keep the maximum 
velocity Vmax)- After Pmax, vehicles can not keep the maximum velocity according to 
their initial gap and the flow begins to decrease although at a slower pace than that 
corresponding to a random initialization. 

Since gapini = L/N — s, the critical density pmax where the transition between 
the free flow and congested state occurs can be calculated as: 

L 

Pmax ~ ' 

By considering, Ax — 2.5 m. s — 2 and setting Vmax — 12 cell/s (this equal 
108 km/k), 

Pmax — 28.57 veh/ km. The associated maximum flow is Qmax — 
Pmax • Vmax = 3085 veh/km. These values are in according to the simulation results 
shown in figure [71 

4. Summary and conclusions 

In this paper a modification of the NaSch model to better capture driver reactions 
to traffic that are intended to preserve safety on the highway was introduced and 
investigated. As a result, three distances that represent the safety distance that a 
driver must have with respect to preceding vehicle if it is going to decelerate, keep 
its velocity or accelerate, were included in the new model. The addition of these 
distances allowed to determine the most appropriate action for a vehicle to undertake 
based on the distance from the vehicle ahead of it and the velocities of the two neighbor 
vehicles. In this way, the velocity setting rules depend not only on the relative velocity 
of neighbor vehicles, they now take into account their relative distance. With this 
modification the normal drivers behavior that bases their driving decisions on both 
relative velocity and relative distance was incorporated. The determination of safe 
distances is in accordance with normal drivers spacing policies, for example, those 
based on vehicle following, constant time headway, etc. Moreover, acceleration and 
deceleration magnitudes for basic maneuvers arc used here arc much closer to those 
reported for normal driving. It is important to emphasize that the order for applying 
rules Rl to i?4 of the proposed model, is not relevant as conditions for each one form 
disjoint sets. Besides, in the model here proposed randomization is applied only to 
vehicles that are in cruising conditions and do not require to brake. In this way double 
braking is avoided, this is consistent with the idea of treating random braking as a 
disturbance that should not occur very often. 

The new model was tested by extended simulations on a one-lane highway with 
a ring topology. The effect of several critical model parameters was analyzed: cell 
lengths, random braking level and initial conditions with the intention of finding 
values that would lead to a proper representation or real traffic behavior. Simulation 
results proved that this new model properly reproduced real traffic flow phenomena. 
Results obtained for homogeneous drivers and different cell lengths exhibit the three 
states observed in real traffic flow: Free-flow states, synchronized states, and stop- 
and-go states, where small jams are present. Synchronized states are associated to 
intermediate-levels of density where vehicles move with a velocity lower than the 
corresponding to free-flow and all vehicles are moving. In this way, stopped cars 
with zero velocity are only present in high density ranges. 

On the other hand, simulation results for R = 0.15 and Ax = 2.5 to, exhibit 
groups of two or more neighbor vehicles traveling at the same speed with the space 
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between them dictated by rules Rl-RA (platoons). For intermediate densities, where 
the synchronized states exist, the number of vehicles that are traveling in platoons 
is about 40%. The platoon formation observed plays an important role in both 
automated and non-automated highway systems. Therefore, results obtained can help 
to elucidate the effects of anticipation coded in the safe distances and the length of 
cell. Besides, the influence of the variation of speed on the flow is also found to be a 
factor of great importance in traffic synchronization. 

Summarizing, the results and discussions in this paper showed the flexibility 
of the CA approach to more complex traffic flow problems. A simple and natural 
modification of the rules of the NaSch model to better capture driver reactions allowed 
to describe the three states of flow observed in the reality. Moreover, the model still 
preserved simplicity of CA model and at the same time, as vehicles' behavior was 
based on a safety analysis to determine the most appropriate action for a vehicle 
to take, made results closer to real highway behavior. Besides, there is an intuitive 
explanation for all rules in the model. It is important to remark that although in this 
paper the model is simulated in a single-lane on a ring, it is possible to apply it to 
complex highway topologies in a satisfactory way. Work in this direction is currently 
in progress. 
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